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We establish the existence of a duality transformation for generic models of interacting fermions 
with two-body interactions. The eigenstates at weak and strong interaction U possess similar 
statistical properties when expressed in the [7 = and U = oo eigenstates bases respectively. 
This implies the existence of a duality point Ud where the eigenstates have the same spreading 
in both bases. Ud is surrounded by an interval of finite width which is characterized by a non 
Lorentzian spreading of the strength function in both bases. Scaling arguments predict the survival 
of this intermediate regime as the number of particles is increased. 
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In noninteracting systems, many-body fermionic states 
are totally antisymmetrized products of one-body states, 
so-called Slater determinants. For strongly interacting 
fermions on the other hand, the eigenstates expressed 
in the basis of Slater determinants have a large num- 
ber of nonzero components, since then the one-orbital 
occupation is in general no longer a good quantum num- 
ber. As the interaction is turned on, the number of such 
components increases, and this crossover from the weak 
to the strong interaction regime has been investigated 
both from the point of view of the statistical properties 
of the spectrum and the eigenstates 0, ||, ^, ^, |[ |j. 
This crossover determines the threshold in excitation en- 
ergy above which Dyson's random matrix theory applies 
for atomic levels (as e.g. in the Ce atom) or nuclear 
energy levels (as obtained in neutron scattering exper- 
iments) and below which, for example, pine quasi- 
particle excitations in chaotic quantum dots are well de- 
fined §]. 

In this article we investigate the reversed crossover 
from strong to weak interaction. Relying both on an- 
alytical arguments and exact numerical diagonalization 
of systems with few particles, we show that the spread- 
ing of the eigenstates over the strong interaction basis 
as the interaction is decreased, is similar to the spread- 
ing over the weak noninteracting basis as the interaction 
is increased. This results in a duality transformation 
between the strong and weak interaction limits. The 
fixed (dual) point of this transformation lies inside a 
finite-width intermediate regime, which is characterized 
by a non-Lorentzian spreading and a maximal complexity 
(which we define below through the structural entropy) of 
the eigenstates in both bases. Scaling arguments predict 
the survival of this intermediate regime as the number 
of particles increases. In this regime, eigenstates proper- 
ties cannot be extracted perturbatively from any of the 
asymptotic limits, but are well captured by our method. 
The eigenstate structure suggests the existence of quasi- 



particle excitations in the strong interaction limit, which 
we were, however, unable to determine. 

We consider the deformed Two-Body Random Ensem- 
ble (TBRE) for n interacting spinless fermions |ll|, [l^ 
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The m different one-body energies are distributed as G 
[—m/2;m/2] so as to fix the mean level spacing A = 1. 
The interaction matrix elements Q^'^ are independent 
random variables with a zero-centered Gaussian distri- 
bution of unit variance. The parameter A = A/{U + A) 
has been introduced to keep the density of states roughly 
constant as U varies. In the U = eigenstate basis, the 
Hamiltonian is represented by a x matrix of size A^ — 
(™) with K = l + n(rn — ri) + ri{ri—l){rn — n)(jn — n—l)/'i 
nonzero matrix elements per row. For a sufficiently large 
number of particles, the many-body Density of States 
(DOS) is well approximated by a Gaussian of width 
Bn - VkA Below we will investigate the prop- 
erties of levels in the middle of the DOS. 

In the U = basis (superscripts (0) and (oo) indi- 
cate the corresponding basis) the eigenstates structure is 
conveniently described by the strength function (SF) 

p(")(i?) = {J2 \i^AiI)\'S{E + Ej- Ea)),, (2) 

A 

and the inverse moments 

€r=((Ei^-4(/)p^)-')^, (3) 



where the averages {...)i,a are taken either over the eigen- 
states (with eigenvalues Ej) of Hq or the eigenstates 
ipA (with eigenvalues Ea) of H. The moment p = 2 is the 
participation ratio (PR) and gives the typical number of 
nonzero components tPa{I)- For model (1) one has ^p^"* S 
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FIG. 1: Width of the strength function p'""' as a function 
of [//A for n e [3,6] and m G [10, 14]. The dotted hue in- 
dicates the Golden rule predicted behavior (^). Upper inset: 
crossover of p'°°^ from delta peak to Lorentzian to Gaussian 
shape for m = 12, n = 5 and U/A = 0.55 (□), 1.09 (+), 
2.15 (x), 4.26 (*), 8.43 (A), 33 (O), and 129.15 (V). The 
solid and dashed lines give a Lorentzian (solid) and a Gaus- 
sian (dashed) fit respectively. Lower inset: strength functions 
p(°°' (□) and p(°) (+) at the dual point Ud « 1.11 for m = 12 
and n — 5. Note that the shape is neither Lorentzian (solid 
line) nor Gaussian (dashed line). 

= NP-^/{2p- 1)!!], where the maximal value 
corresponding to Dyson's Gaussian Orthogonal Ensem- 
ble (GOE) is achieved for U = oo [^. There are, beside 
the one-body spacing A, two important energy scales [Q] : 
the mean spacing between states directly coupled by the 
two-body interaction A^"'' — B^^ /K 4A jmn} and the 
n-body spacing A^"-' = Bn'^ /N , where Bn^ « n(m — 7i)A 
is the n-body band (and hence B^^ « 2(to — 2)A). As U 
increases, it has been found that p'^'^^E) undergoes two 
transitions from a delta peak to a Lorentzian shape first, 
then to a Gaussian shape [|l], |l2[ |lj]. In the Lorentzian 
regime, the width r(*" of p(")(£:) is given by the golden 
rule r("' cx C/VAi*" and the PR is ~ r(o)/Al°^ §. 
The Lorentzian regime is defined by the two conditions 
> 1 and r(o) < B'i? (l|, 0. In the dilute limit 
1 ^ n ^ TO, these conditions translate into 

A/V^ <t:u < A/y^. (4) 

When the Gaussian regime is entered, the SF spreads 
over the full bandwidth so that the eigenstates have 
a finite fraction of nonzero components ^j""* — 0{N). 
Having summarized some of the known results for the 
eigenstates structure as the interaction is switched on, 
we now turn our attention to the reversed problem and 
decrease U starting from U — oo. Both the SF 
and the moments can be defined in the same way 
as above provided the U = oo basis is fixed individually 
for each realization of the Hamiltonian (1). Since the 



occupation operator = c\ca does not commute with 
-ffi, Hq induces one-body transitions between different 
eigenstates of Hi , leading to an increase of ^1°°^ > 1 and 
a broadening of as U decreases. The ilo-induced 
transitions also lead to two successive crossovers of the 
SF p^°°\ first from a delta peak to a Lorentzian shape, 
then to a Gaussian shape, and this is shown on the up- 
per inset to Fig. ^ As is the case for r(°), the Golden 
rule gives a good estimate of the width r(°°) of the SF 
in the Lorentzian regime, as we now proceed to show. 
Under the assumption that for [/ = oo, the TBRE has 
random ergodic eigenstates in the middle of its spectrum 
JlSf , all eigenstates are directly connected to each other 
by Ho so that Ai°°^ = Ai°°^ = Bi°°^/N. In the Hmit 
of large number of particles and orbitals, the U = oo 
density of states is well approximated by a Gaussian 
with a width bIT^ « VkAU « n(m - n)A 0, ||, 
and one has A^°°'' « An"-* cx n^m — n)A/N . Assuming 
random ergodic many-body eigenfunctions at ?7 = oo 
I pd] , , the transitions matrix elements have a variance 
{Hi) = Y.o.{^o./UflN « (n(TO - n)A-'/UflN. The 
Golden rule then predicts 

r(°°) ocn(m-n)AVC/^. (5) 

This prediction is confirmed by numerical data presented 
in Fig. 1^. These data, as well as those to be presented 
below, have been obtained via exact diagonalization of 
systems of up to = 3432 (corresponding to n = 7 and 
TO — 14), performing averages over 20 (for the largest 
systems) to 500 realizations of Hamiltonian (1) for each 
parameter set and avoiding the tails of the DOS by keep- 
ing only 33% of the states in the middle of the spectrum. 

In the Golden rule regime, the PR is given by S^^^ = 
r(°°V^"°°^ ~ N{A/Uy and according to the conditions 
^!^°°^ > 1 and r(°°) < Bi°°\ the Lorentzian regime in the 
U = CO basis is bounded by the inequalities 

A<U < VaA. (6) 

This is confirmed qualitatively by the upper inset to 
Fig. |l| and by numerical data to be published elsewhere 

The remarkable fact that Ai°"^ and An"^ have the same 
parametric dependence in n and m implies that P'"' = 
r(°°) and ^!^^ = ^2°°^ can be both satisfied for Ua - 
Ajv}!^. This is illustrated on the lower inset to Fig. |l| 
and the inset to Fig. ||. Remarkably enough, we have 
found that at Ud^ the PR's take a universal value — 
^2°°-' sa 0.8 ^2"*^^ independently on n and m. This calls 
for a rescahng \J U/Ud for and U -> Ud/U 

for ^^°°VC2'°^ after which all data for the PR faU on top 
of each other, as can be seen on Fig. ^ 

The inset to Fig. ^ shows that at Ud, the higher mo- 
ments with p > 2 also cross. This behavior has been 
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FIG. 2: e2°VC°°'' (empty symbols) and ^^°°VC2'°'' (Ml 
symbols) vs. the rescaled interactions U/Ud and Ud/U re- 
spectively. Different symbols correspond to different n G [2, 7] 
and m £ [10, f 6]. Inset: fZ-dependence of 5p°' (empty sym- 
bols) and Cp°°' (full symbols) for p = 2 (□), 3 (Q), 4 (A) and 
5 (V), for m = 12 and n = 5. The dashed line indicates the 
(p- independent) intersection point Ud/^ ~ 1.11. 



found to hold for all m and n. It implies that the eigen- 
states have exactly the same spreading over both bases, 
and that the SF's have not only the same width, but 
also exactly the same shape at Ud (^~^ are moments of 
the SF), as shown in the lower inset to Fig.|l]. More- 
over, we found that at the dual point the SF satisfies 
the scaling T{Ud) oc mn'^/^A as follows from (|^) and 
Ud (X A/n^/". We then calculate the PR ^^''^ in the 
U ~ Ud (dual) basis, and numerical results are shown in 
Fig. ^. The symmetry around U ^ Ud is evident, and for 
both [/ — !■ and t7 ^ oo limits, there is saturation at 
0.8 '^2^'^^. These results establish the existence of a 



Ad) 

duality transformation U* = U'^/U connecting the strong 
and weak interaction regimes, so that the PR's and 



the SF satisfy S^^^U) = Cr\U 
p(o){U) ^ p^°^HU*) and p('^)({7) = 



f(oo) , 



p('^Hu*). 



The above estimate for Ud ^ A/n^/^ relies on the 
assumptions that the Golden rule correctly estimates 
the width of the SF in both bases, and that = 
p(o,oo)/^(o,oo)^ Equivalently, this requires to be in the 
Lorentzian regime, which is however impossible as the 
two conditions (^) and (p|) are mutually exclusive, and 
the predicted Ud ~ A/ti^'"* lies outside both Lorentzian 
regimes. One may thus wonder if its predicted paramet- 
ric dependence in m and n makes any sense at all. Ud 
as a function of both m and n is shown in the insets 
to Fig. ^. First, it is seen that the m-dependence of Ud 
is very weak (we found Ud ~ Ato~" with a < 0.08 for 
a range m S [8,20]). Even though quite weak, this m- 
dependence presumably indicates that we are not deep 
enough in the dilute limit. Second we extracted the n- 
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FIG. 3: C^''VC?°^ calculated in the U ^ Ud basis for m = 12 
and 71 = 5. Left inset: n-dependence of Ud for m — 14. Right 
inset: m-dependence of Ud for n — 3 (O) a-nd 4 (□). 
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FIG. 4: Structural entropy Sstr = S - ln(^2) in the !7 = 
(empty symbols) and U = oo (full symbols) basis for m = 14, 
andn = 2(n),3(0),4(A), 5(V). 



dependence of Ud for m = 14 where we had the largest 
range n S [2,7]. We got Ud ^ An~^, with an exponent 
P e [0.3; 0.5]. These bounds on (3 are compatible with 
the Golden rule estimates which exclude Ud from both 
Lorentzian regimes, i.e. 



A/V^ < C/d < A, 



(7) 



giving /3 e [0,0.5]. We conclude that Ud ^ A/n^/'* is 
in good qualitative agreement with our numerical re- 
sults, and that the inequalities (^ define an intermediate 
regime which is nonperturbative in the two bases, and 
whose width increases parametrically with n. 

Additional structures in the wave functions can be 
captured by the structural entropy Sgtr, which is de- 
fined by subtracting from the Shannon entropy S 



(E,I^A(/)pln|V'A(/)P>A 



the log of the PR: 3^°;°°^ 
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5(0,00) _ in^(o,oo) ^ It measures the contribution to 
the Shannon entropy which is not contained in the PR, 
and thus not in the bulk of the SF. Numerical results are 
presented in Fig. ^. Expectedly, the crossing of S'f°i.'°°'' 
occurs at Ud (since all the moments determine Sstr) 
but the crucial feature of Fig. ^ is that almost imme- 
diately after the crossing, S'str takes on its asymptotic 
GOE value 5,^°'^ = 0.3689 ... - and therefore Ua is the 
crossover point to the regime with GOE-like behavior 
in both bases. The second crucial feature of Fig. |4| are 
the peaks in S'str located on both sides of Ud- Such peaks 
have already been observed in various models where they 
indicate unusually large fluctuations of eigenstates p8| . 
It is expectable that these peaks extend over the interme- 
diate, non-Lorentzian regime AZ-^n < U < A, which is 
borne out by the data plotted in Fig. 4, where once S'str 
is plotted against U/Ud — Un^^'^/A, the peaks broaden 
symmetrically as more particles are added. 

Examples of models where such a duality between lo- 
calized (with ^2 = 1) and delocalized (with ^2 = 0{N)) 
asymptotic regimes is related to a sharp metal-insulator 
transitions are provided by one-dimensional lattices with 
quasi-periodic potential (where the duality connects mo- 
mentum and spatial eigenfunction coordinates) [^9| and 
tight-binding models for strongly correlated fermions 
pO| . Our model is however fundamentally different in 
that the duality point Ud is protected by a finite-sized 
interval characterized by a non-Lorentzian spreading of 
the eigenstates over a finite fraction of both the strong 
and weak interaction eigenstate. According to (|^), this 
intermediate regime survives as m and n increase and this 
results in a smooth crossover (and not a sharp transition) 
to fully developed chaos, characterized by the collapse of 
all the data points over one single curve shown on Fig. ^ 
(and not a single crossing of the curves at U/Ud = 1). 
This is to be put in perspective with the recent contro- 
versy over whether or not particle-hole excitations un- 
dergo a delocalization transition in Fock space similar 
to the Anderson transition as their excitation energy in- 
creases ^ |9| (increasing the excitation energy reduces 
the level spacing and is thus equivalent to increasing the 
interaction strength U). The results we presented cor- 
roborate the conclusions drawn in [Q of the existence of 
a smooth crossover and not a sharp transition. 

In summary we have established the existence of a du- 
ality transformation between the weak and strong inter- 
action regimes of deformed TBRE. At the duality point 
Ud, the eigenstates have GOE fluctuations (for U > 
{<)Ud in the U = (00) basis) and Ud is surrounded 
by a finite-sized, intermediate regime where both limit- 
ing bases fail to describe the eigenstates. Together with 
the Lorentzian form of p'^^^ at large U ^ A, this du- 
ality suggests the existence of quasiparticle excitations 
at [/(J ^ [/ < cxD, which we were however not able to 
identify. The existence of this duality transformation is 
related to the two-body nature of the interaction. For 



fc-body interaction, n ^ k > 2, the U = 00 bandwidth 
is given by Bi°°^ - An'^/^m''/^ > si"^ |ll|, |l|, and 
the PR in both bases are equal to each other for Ud oc 
Ami-'=/2n3/4-fe/2. We then have T^°\Ud) cx Amn^/^ 
while r(°°)(f7d) (X Am'=/2ni/2+fe/2 go that the SF's should 
differ at the intersection point of the PR's for fc 7^ 2. 
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